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ABSTRACT 

We have investigated the effect of spatial resolution on determining pencil-beam like 
velocity widths and column densities in galaxies. Three 21-cni datasets are used, the 
HIPASS galaxy catalogue, a subset of HIPASS galaxies with ATCA maps and a high- 
resolution image of the LMC. Velocity widths measured from 21-cm emission in local 
galaxies are compared with those measured in intermediate redshift Damped Lyman-a 
(DLA) absorbers. We conclude that spatial resolution has a severe effect on measur- 
ing pencil-beam like velocity widths in galaxies. Spatial smoothing by a factor of 240 
is shown to increase the median velocity width by a factor of two. Thus any differ- 
ence between velocity widths measured from global profiles or low spatial resolution 
21-cm maps at z = and DLAs at z > 1 cannot unambiguously be attributed to 
galaxy evolution. The effect on column density measurements is less severe and the 
values of cLNdlaI dz from local low-resolution 21-cm measurements are expected to 
be overestimated by only '--^lO per cent. 
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1 INTRODUCTION 

Damped Lyman-a (DLA) absorption systems provide the 
, most practical measure of neutral gas in intermediate red- 
shift galaxies. However, measurements of DLAs provide in- 
formation in one dimension only - the redshift axis, and 
only along a single pencil-beam line-of-sight. In rare cases 
of close binary QSO pairs or gravitationally lensed QSOs, 
spatial information on the absorbers can be deduced (e.g. 
iRauch et ai]l2002l : lEUison et al.ll2004 . The nature of DLAs 
can be inferred by comparing observations of DLA velocity 
widths to models of galaxy evolution. 

DLA kinema tics and models of galaxy evolution were 
first compared bv lProchaska fc Wolfd ^ll997^ who concluded 
that rapidly rotating protogalactic discs provide the best fit 
to DLA kinematics. If DLAs are the progenitors of spiral 
discs, then their velocity widths are expected to match that 
of present day galaxies. This interpretation however seems 
to run contrary to cold dark matter theory. 

The current structure formation paradigm predicts that 
star formation first occurs in subgalactic clumps housed 
in the most massive dark matter halos. This framework 
is well matched by the assembly of stars and dark mat- 
ter into galaxies, where more massive halos transform a 
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larger fraction of t heir gas into stars at high redshift 
i Jimenez et aljl2005l) . Further evidence is provided by the 
clustering strengths of Lyman break galaxies (LBGs) at 
2 ~ 3, showin g that LBGs are embe dded in massive dark 
matter halos jAdelberger et all 12005*). DLAs on the other 
hand are quite distinct from LBGs (Hopkins ot al. 2003), 
perhaps corresponding to lower-mass systems with lower ro- 
tational velocities. Since lower-mass systems convert their 
gas into stars at a slower rate they contribute significantly 
to the total DLA cross-section (averaged over time). 

Hydrodynamic simulations demonstrate that irregu- 
lar protogalactic clumps can equally reproduce the ob- 
served velocity width distribution of DLAs through a com- 
bination of rotation, random motions, infall and merging 
iHaehnelt et alJll998l) . Both semi-analytic and numerical 
models show that the average velocity width o f DLAs is 
expected to increase with evolution to g = iCen et alJ 
l200^:'Okoshi et al.'2004': lNagamine et all2004) . Other work 
(Mai ler et al. 2001) interprets the kinematics of DLAs as ev- 
idence that they arise in galaxies undergoing tidal stripping 
and mergers. 

Because very few DLAs are known at low redshift, it 
is difficult to measure their statistics directly. Fortunately, 
at 2: = the alternative method of 21-cm emission can be 
used. It has the added advantage of known galaxy properties 
so the distribution of DLA galaxy properties such as lumi- 
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nosity and H I mass can be inferred. iRosenberg fc SchneideJ 
i2003r have compared global 21-cm velocity widths to ve- 
locity widths of intermediate redshift DLAs. They conclude 
that evolution is not the likely cause of the difference be- 
tween the two populations since velocity widths in their 
sample of local galaxies are as large or larger than the DLA 
velocity widths. Instead they suggest that the difference be- 
tween the two populations could be caused by small number 
statistics or differences in the way the two populations were 
measured. Unfortunately 21-cm emission measurements al- 
ways come with the caveat of spatial resolution as the typ- 
ical beam size is much larger than the region probed by 
DLA observations. The extent of the absorbing region in 
DLAs is only as large as the projected angular extent of the 
UV-emission region of the background QSO. An upper limit 
to the optical emission region of one QSO is known to be 
<0.03 pc by microlensing ( Wyithe et al. 2002). The effec- 
tive angular extent of absorption in each DLAs is expected 
to be roughly constant with redshift as the angular diameter 
distance varies less than 20 per cent over 1 < z < 4, where 
most QSO probes and DLAs are detected. 

In this paper the relationship between H I velocity 
widths in local galaxies and intermediate redshift DLAs is 
investigated. Does the velocity width distribution of neutral 
gas in galaxies change with evolution to 2: = ? In address- 
ing this question, the relationship of pencil-beam DLA mea- 
surements to the global H I profile of galaxies and the effect 
of spatial resolution on determining H I column densities and 
velocity widths in 21-cm maps must also be investigated. 



2 DLA KINEMATICS 

Here we test the evolution of galaxy kinematics by com- 
paring DLAs at intermediate redshifts with gas-rich galax- 
ies at the present epoch. Velocity widths from two different 
21-cm datasets are utilised: the 4315 galaxies in HICAT, 
which is t he H I Parkes All Sky Survey (HIPASS) galaxy 
catalogue iMever et al.ll2004) and a subset of 34 HICAT 
galaxi es with Australia Telesco pe Compact Array (ATCA) 
maps iRvan- Weber et al.ll2003() . Velocity widths from these 
21-cm observations are comp ared with a c omplilation of ve- 
locity widths from 94 DLAs JWolfe et alj EoOS^. The veloc- 
ity widths of DLAs are measured using low ionisation metal 
lines such as Si^, Fe"*" and Ni^ to trace the neutral g 
the Lya line, by definition, is saturated and does not re- 
solve the DLA kinematics. Velocity profiles of low ionisation 
metal lines trace each other well and therefore are assumed 
to trace the kinematics of the cold gas. This assumption 
could be v erified with high re solution lUE spectra towards 
SN 1987A iBlades et alJll988h . 

The global Hi profile statistic used from HICAT is 
W5Q, the velocity width at 50 per cent of the profile max- 
imum. Comparisions are made with W50/2 as a single sight- 
line though a galaxy will presumably only sample one 
half of a galaxy's rotation. The relative contribution of 
each galaxy has been weighted by the H I mass dependent 
DLA cross-section (.Rvan-Wc bcr et al. 2003) . Figure^shows 
the distribution of W50/2 and velocity widths of low ion- 
isation lines in DLAs. A Kolmogorov-Smirnov (K-S) test 
shows the probability that the two data sets are drawn 
from the same distribution is Pks ~ 0.01, compared with 
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Figure 1. Velocity width histogram comparing w)5o/2 from HI- 
CAT galaxies to low ionisation metal lines in DLAs. 

IRosenberg fc Schneideil (l2003h . who find Pks = 0.89 and 
conclude that evolution in the kinematics of H l-rich galax- 
ies cannot be ruled out (a value of Pks ^ 0.99 would rule 
out evolution). The same conclusion is then reached with 
the HICAT dataset: evolution cannot be ruled out, the two 
populations may differ due to evolution, but could also difi'er 
due to small number statistics or differences in the way the 
velocity widths are measured. 

Differences in the way the velocity widths are mea- 
sured can be reduced by measuring the 21-cm widths with 
a method that most closely resembles the determination 
of QS O absorption line velocity widths. IProchaska fc Wol^ 
il997ll measure the velocity width across the central 90 per 
cent of the optical depth. The ATCA 21-cm maps are used 
to measure a pixel velocity width along the velocity axis of 
each spatial pixel in each galaxy map. Only spatial pixels 
that satisfy the DLA cutoff, A^m^ 2 x 10^° cm~^, are in- 
cluded. The velocity width is determined for the central 90 
per cent of the integrated profile, this measurement is de- 
noted by AvgQ%. The average restored beam of the ATCA 
maps is 1', corresponding to a physical extent of 2 to 7 kpc. 

Figure |5| shows that neither W50 nor U150/2 provides a 
good statistical representation of the median Awgo% velocity 
width. Only 18 of the 34 galaxies in the sample have W50/2 
measurements that fall within 30 km s~^ of their median 
Aw9o% . This is not a result of the DLA nor 90 per cent cutoff, 
as a similar spread of median velocity widths is found when 
all spatial pixels with Nm^ 10^^ cm^^ or 100 per cent of 
the integrated profile is considered. 

The reason why w^o widths are not well represented 
by pixel velocity widths is due to the variety of galaxy 
types in the sample. This is problematic as it is now well 
established that the local galaxy population weighted by 
DLA Hi cross section s ar nples a var iety of galaxy types 
iRvan- Weber et al.l2003l : IZwaan et al]l2005l) . Low mass and 
dwarf galaxies tend to have Gaussian velocity profiles that 
remain approximately constant over the area of the 21-cm 
map. The AvgQ% method acting on a Gaussian profile, will 
return velocity widths which are by definition 29 per cent 
larger than the FWHM, that is, ~ 60 per cent larger than 
W50/2. Face-on spiral galaxies are expected to follow a sim- 
ilar overestimate of velocity widths, but not as extreme as 
the low mass galaxies. On the other hand, the pixel veloc- 
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Figure 2. Global w^(j/2 velocity width versus median Avqq% for 
ATCA galaxies. The dotted line is where the two axes are equal. 
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Figure 3. Histogram comparing At)go% from the ATCA galaxy 
sample with DLA velocity widths. 



ity widths of large edge-on spiral galaxies that have many 
pixels spread over a large area are more likely to be lower 
than u)5o/2. These trends are summarised in Figure|5| As ex- 
pected, all the low mass (log Mm/ Mq< 9) galaxies lie on or 
above the line where W50/2 equals the median /S.Vqq%. Most 
face-on (inclination < 45°) high mass (log M-ai/MQ> 9) 
galaxies also lie on the same side of the line. However, most 
massive edge-on (inclination > 75°) galaxies have median 
Ad9o% that are well below their global W50 /2 measurements. 
The range of galaxy types means than any statistic of the 
pixel velocity width, e.g. some fraction of the maximum ve- 
locity width, will not be well represented by any global veloc- 
ity width statistic. We conclude that the global 111^,0 value 
does not provide a good estimate of pencil-beam velocity 
widths, as measured in DLAs. 

Pixel velocity widths from the ATCA galaxies are com- 
pared with those from DLAs in Figure O The K-S proba- 
bility is calculated to be less than 0.001, indicating that the 
two distributions are different. Does this mean that there is 
genuine evolution in velocity width from the intermediate 
redshift DLAs to present day H l-rich galaxies, or is the dif- 
ference still due to the measurement techniques? The lack of 
high velocity widths could be caused by treating galaxies as 
individual systems, if their satellites were considered, larger 



velocity widths could be invoked. A deficit of low velocity 
widths cannot be caused by velocity resolution since the 
resolution is 3.3 km s~^ in the ATCA images. The galaxies 
that tend to have the lowest velocity widths also have the 
smallest contribution to dNoLA/dz, however there seems to 
be a complete lack of Avr,o% values below 40 km s~^. The 
DLA column density cutoff is well above the noise limit in 
the ATCA images, so the signal-to-noise ratio cannot be an 
issue. Discrepancies between the velocity width of gas mea- 
sured in emission and absorption would ideally be tested 
with UV spectroscopy towards QSOs rece ntly id entified be- 
hind the Large Magel lanic Cloud (LMC, ICeha e t al. 200l 
iDobrzvcki et alj|2005l) . an experiment that awaits the next 
space-based UV spectrograph. 



3 VELOCITY WIDTH STATISTICS WITH 
THE LMC 

The absence of low velocity line widths could be caused by 
the relatively low spatial resolu tion of the AT CA galaxy 
maps. Following on from iProchaska et al.l ll2002ll a compari- 
son can be made using the H I data cube of LMC. The LMC 
data cube is a high-resolution 21-cm data cube con structed 
from t he combined Parkes and ATCA survey of Ki m et alJ 
i2003(l . The data has a velocity resolution of 1.6 km s~^ and 
a column density sensitivity of 4 x 10^^ cm~^ (for a line width 
of 40 km s~^). The LMC has an Hi mass of 3.8x10*^1/0, a 
global Hi velocity width, 1050=80 km s^^, is located at a 
distance of 50 kpc (Stavclcv-Smith ct al.^20Q3) and has an 
incfination of 42° IWeinberg fc Nikolae^llioOll) . The ATCA 
galaxy and LMC maps both have an angular resolution of 
~l'. However, the LMC is 80 to 480 times closer than the 
galaxies in the ATCA sample, leading to a large difference 
in physical spatial resolution, 2 — 7 kpc per beam for the 
ATCA galaxy sample compared with 15 pc per beam for 
the LMC. 

Since the LMC data cube contains over 4 x 10® spatial 
pixels, a random sample of 10* pixels is used in this analysis. 
The Aw9o% method was applied to this random sample of 10* 
pixels and the results are given in Figure|l| Only pixels that 
satisfy the DLA cutoff, A^hi ^ 2 x 10^"cm~^, were included. 
The LMC cube was convolved with a 2-D circular Gaussian 
in the spatial plane. The LMC cube was smoothed on three 
scales such that the resolution in the smoothed cubes were 
80, 240 and 480 times greater than the original. A resolution 
factor of 240 corresponds to a final resolution of 4°, effec- 
tively placing the LMC at 12 Mpc, or the maximum distance 
at which galaxies with Hi masses similar to the LMC are 
found in the ATCA galaxy sample. 

The velocity width histogram for this smoothed cube is 
plotted in Figure |1] The difference is stark and shows that 
the median velocity width from the smoothed cube using 
the A«go% method is now comparable to wso, rather than 
1050/2. Decreasing the resolution shows an overall increase 
in the median velocity width is because the velocity width 
of the smoothed set of pixels will return the maximum ve- 
locity width rather than the average velocity width. The 
median velocity width continually increases with decreas- 
ing resolution: the unsmoothed cube with a resolution of 1' 
and smoothed cubes with resolutions 1.3°, 4° and 8° have 
median velocity widths of 41, 64, 81 and 87 km s~^ respec- 
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Figure 4. Histogram of AvgQ% produced from 10'* random pixels 
in the LMC cube satisfying A^hi > 2 x 10^° cm"^. The red his- 
togram is the same calculation for the LMC cube convolved with 
a circular Gaussian to produce a smoothed image with FWHM = 
4° . The dotted and dashed lines are U150 and W50/2 respectively. 



Figure 5. Histogram of log column density from 10'' random pix- 
els in the LMC cube. The red histogram is of the same smoothed 
image described in Figure]^ No A^jji restrictions are applied to 
either data set. 



lively. Turning this statement around, as the resolution in- 
creases, the distribution of velocity width spreads from half 
the low resolution minimum velocity width to the low reso- 
lution maximum velocity width and the median is reduced 
by half. Applying this general trend to the ATCA galaxy 
sample in Figure |3] a distribution more like the DLAs is ob- 
tained. Therefore this resolution effect cannot be decoupled 
from evolution in the velocity width distribution between 
local gas-rich galaxies and higher redshift DLAs. 

The highest resolution LMC cube probes structure on 
scales 500 times greater than that of DLA measurements. 
Decreasing the spatial resolution in the LMC cube a further 
240 times essentially doubles the median velocity width. If 
Hi structure is self-similar on smaller scales, then one would 
expect the velocity widths to decrease further as the scale 
on which DLAs are measured is approached. 



4 THE EFFECT OF SPATIAL RESOLUTION 
ON COLUMN DENSITY 

Since spatial resolution was found to have such a dramatic 
effect on velocity width measurements, what is its effect on 
the other 21-cm based quantities? There is only one case of 
H I emission detected from a DLA galaxy. The DLA galaxy 
SBS 1543+593 (Bowcn ct al. 2001) has been mapped using 
GMRT IChengalur fc Kanekaill2002i) in 21-cm. The column 
density measured at the QSO location in the GMRT map is 

5 X lO^^cm"^, in reasonable agreement with the Lya column 
density of 2.2 x lO^'^cm"^. Here the effect of resolution on 
the column density distribution function, f{Nm), and the 
cross-section of Hi satisfying the DLA criteria, that is the 
number density of DLAs per unit redshift, dNoLA/dz, is 
investigated. 

The column density histogram for 10"* random spatial 
pixels in the LMC cube is given in Figure|H|together with the 
histogram derived from the LMC cube convolved with a cir- 
cular Gaussian to produce a cube smoothed in the spatial 
plane with FWHM — 4°. The full range of column densi- 
ties are considered here. Lowering the resolution in the cube 
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Figure 6. LMC fiNm) (black squares), with arbitrary normali- 
sation. Green circles and red triangles give the LMC /(A'^hi) for 
smoothed cubes with FWHM = 1.3° and 4° respectively. The dot- 
ted line is the double power law fit from Ryan- Weber et al. (2003), 
with normalisation correction from Ryan- Weber et al. (2005). The 
solid blue line is /{Nm) from Zwaan et al. (2005). 

means that any fine structure in low or high column density 
regions is washed out. The histogram refiects this qualita- 
tive expectation and shows that the distribution of column 
densities are truncated at both the low and high ends. 

The column density distribution function will naturally 
be influenced by this change in di s tribut ion. F igure El shows 
fJNm) from iRvan- Weber et alJ (|20o3) and IZwaan et all 
||2Q23) together with f{Nm) calculated from 10* random 
spatial pixels in the LMC cube at the full resolution and 
smoothed at two different levels. Since only one galaxy, the 
LMC, is considered, the same arbitrary normalisation is 
applied to aU the LMC f{Nm) points. None of the LMC 
f{Nm) functions resemble the shape of f{Nm) calculated 
from many galaxies. The shape of f{Nm) does vary be- 
tween galaxies and a turnover at A^hi~ lO^'^ cm"'^ occurs 
in ~ 10 per cent of the ATCA gala:xies (see Ryan- Weber et 
al. (2003) figure A2). The plot shows that decreasing the 
resolution creates a hump in the centre of the function. Ap- 
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plying this trend to f{Nm) calculated from the lower resolu- 
tion ATCA maps, suggests that the 'real' underlying f{Nm) 
could be described by a single power law, as the lower A^'hi 
slope would steepen and the upper A^hi slope would flatten. 
Although this change in shape is based on one galaxy only. 

The calculation of the number of DLAs per unit red- 
shift at z — 0, dNoLA/dz, is also affected by resolution. 
The DLA cross-section of the LMC increases slightly with 
decreasing resolution: the increase in the DLA cross-section 
from the original resolution to cubes smoothed to a resolu- 
tion of 1.3, 4 and 8° is 11, 12, and 12 per cent respectively. 
Hence the overall effect of resolution on local dNpLA/dz cal- 
cula tions ( Rosenberg & Schneider 20031: iRvan- Weber et all 
l200a..2005; Zwaan et al...2005,) is small, and this 11-12 per 
cent increase lies within the uncertainties for the measured 
values of dNoLA/dz. 



5 SUMMARY 

The results highlight the difference in conclusions when ve- 
locity widths are measured in different ways. The w^o ve- 
locity width from a galaxy's global If I profile is not well 
represented by the distribution of pixel (pencil-beam like) 
velocity widths in a galaxy. Instead, the relationship be- 
tween wso and median pixel velocity width was found to 
depend on galaxy mass and inclina tion. This work provides 
a solid basis for the suggestion by iRosenberg fc Schneider! 
i2003h that evolution in DLA kinematics (measured by 21- 
cm global velocity widths at z = 0, and absorption line 
profiles at higher redshift) cannot be decoupled from the 
different measurement techniques. Thus using DLA velocity 
widths to infer the total velocity width of a galaxy and re- 
lated quantities, e.g. mass, must also be used with caution 
(e.g. Godfrey & Francis 2005,). 

Furthermore, when we use a method that most closely 
resembles that of DLA velocity width measurements, the 
Aw9o% method, the same conclusion is reached. The distri- 
bution of velocity widths measured in DLAs at intermediate 
redshifts is different from those measured along the line-of- 
sight of individual spatial pixels within galaxies. It is not 
possible to decouple the effect of spatial resolution in low- 
resolution 21-cm maps from evolution in galaxy kinemat- 
ics. Unfortunately the large difference in spatial resolution 
sampled by 21-cm emissi on and Lya abso r ption cannot be 
reconciled. The results of IProchaska et all (|20o3) still hold, 
as a single low-mass galaxy such as the LMC cannot ac- 
count for high (> 100 km s~^) DLA velocity widths. How- 
ever differences still exist between the highest resolution 21- 
cm data, which have multi-component profiles overlaping in 
velocity coverage, and DLA profiles, which tend to decom- 
pose into several components distinct in velocity coverage 
(e.g. Wolfe et al. 2005). 

Column density measurements and related calculations 
are not adversely affected by spatial resolution. Lowering the 
resolution in 21-cm data means that fine structure in low or 
high column density regions is washed out. This causes a 
central hump in the column density distribution function. 
The expected increase due to spatial resolution effects in 
ca lculating dNnLA/dz at z=0 is w ithin uncertainties quoted 
in IRosenberg fc Schneideil J2003D . iRvan- Weber et alJ J2003l 
|2005^ and lZwaan et alJ ^2005^ . 
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